Non-Poisson processes: regression to equilibrium versus equilibrium correlation 

functions 
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We study the response to perturbation of non-Poisson dichotomous fluctuations that generate 
super-diffusion. We adopt the LiouviUe perspective and with it a quantum-hke approach based on 
sphtting the density distribution into a symmetric and an anti-symmetric component. To fit the 
equiUbrium condition behind the stationary correlation function, we study the time evolution of the 
anti-symmetric component, while keeping the symmetric component at equilibrium. For any realistic 
form of perturbed distribution density we expect a breakdown of the Onsager principle, namely, of 
the property that the subsequent regression of the perturbation to equilibrium is identical to the 
corresponding equilibrium correlation function. We find the directions to follow for the calculation 
of higher-order correlation functions, an unsettled problem, which has been addressed in the past 
by means of approximations yielding quite different physical effects. 

PACS numbers: 



I. INTRODUCTION 

Complex physical systems typically have both nonlin- 
ear dynamical and stochastic components, with neither 
one dominating. The response of such systems to ex- 
ternal perturbations determines the measurable charac- 
teristics of phenomena from the beating of the human 
heart to the relaxation of stressed polymers. What dis- 
tinguishes such complex phenomena from processes suc- 
cessfully studied using equilibrium statistical mechanics 
is how these systems internalize and respond to environ- 
mental changes. Consequently, it is of broad interest to 
determine which of the prescriptions from equilibrium 
statistical physics is still applicable to complex dynami- 
cal phenomena and which are not. Herein we address the 
breakdown of one of these fundamental relations, that be- 
ing, the Onsager Principle. In the case of ordinary statis- 
tical mechanics an exhaustive treatment of the relaxation 
of perturbations to equilibrium can be found in Ref. [l| . 
Let us consider as a prototype of ordinary statistical me- 
chanics the case when the stochastic variable under study 
^(t) is described by the linear Langevin equation 

§ = ~im+vit), (1) 

where the random driving force ri(t) is white noise. Let 
us imagine that ^(i) is the velocity of a particle with 
unit mass and a given electrical charge. Furthermore, we 
assume that this system reaches the condition of equi- 
librium, and at a given time i = 0, we apply an elec- 
trical field E{t). The external field E{t) is an arbitrary 



function of time, fitting the condition that E{t) — 0, for 
t < 0. The adoption of linear response theory yields the 
prescription for the mean response of the system to the 
external field 



im) 



<^^{t')E{t~t')dt\ 



(2) 



where $^(t) is the equilibrium correlation function of ^. 
In this familiar case, of dissipative Brownian motion, the 
autocorrelation function of the particle velocity is the 
exponential exp(— 7^). 

There are two limiting cases of time dependence of the 
perturbation: (a) the electric field is proportional to the 
Heaviside step function, E{t) — KQ{t), and is therefore 
a constant field after it is turned on at i = 0; (b) the 
electric field is proportional to the Dirac delta function, 
E{t) — K5{t), and is consequently an initial pulse that 
perturbs the particle velocity. In these two limiting cases 
we obtain for the velocity of the Brownian particle 



{m)=K ^dt')dt' 



and 



(C(t)) = K^iit), 



(3) 



(4) 



respectively. These two limiting cases show that in the 
case of ordinary statistical mechanics the system's re- 
sponse to an external perturbation is expressed in terms 
of the unperturbed autocorrelation function. We shall 



refer to the conditions of Eq. (jJJ and Eq. Q as the 
Green-Kubo relation and the Onsager relation, respec- 
tively. 

The search for a dynamical derivation of anomalous 
diffusion, that being where the mean square value of the 
dynamic variable is not linear in time, has been a sub- 
ject of great interest in recent years. There are two main 
theoretical perspectives on how to explain the origin of 
anomalous diffusion. The first perspective is based on 
the assumption that there are unpredictable events, that 
the occurrence of these events obey non-Poisson statis- 
tics, and is related to the pioneering paper by MontroU 
and Weiss ^. The other perspective rests on the as- 
sumption that the single diffusion trajectories have an 
infinite memory. The prototype of the latter perspective 
is the concept of fractional Brownian motion introduced 
by Mandelbrot y] ■ A problem worthy of investigation is 
as to whether or not, in the case of anomalous diffusion, 
the response to external perturbation departs from the 
predictions of Eqs. Q and 10} . In the last few years, 
this problem has been addressed by some investigators 
0,|l,|a,0lljl3- These authors have discussed the Green- 
Kubo relation of Eq. Q . Notice that in the special case 
of ordinary statistical mechanics this relation can also be 
written in the following form 



{x{t)) 



K 

W) 



{x\t)), 



(5) 



To understand how to derive this equation, originally 
proposed by Bouchaud and George |l3 , we have to refer 
ourselves to the following equation of motion 






(6) 



Since ^(t) is a fluctuating velocity, it generates spatial 
diffusion and we denote by x(t) the position of the cor- 
responding diffusing particle. The external field affects 
the velocity fluctuation and, consequently, the diffusion 
process generated by these fluctuations. In the absence 
of perturbation, the second moment of the diffusing par- 
ticle, {x{t)'^)o, obeys the prescription 

{x{tf)o = 2{e) f dt' f dt"^^{t"). (7) 

Jo Jo 

It is straightforward to prove that Eq. (0) yields Eq. jSJ. 
This is done by considering the mean value of Eq. (|5|) 



d{x) 
dt 



- m))- 



(8) 



The time integration of the left hand term of this equa- 
tion yields the left hand term of Eq. Q and the time 
integration of the right hand term of it, using Eq.Q and 
Eq.Q, yields the right hand term of Eq. ©. 

The relation of Eq. jSJ is denoted as generalized Ein- 
stein relation^ because it might hold true also when the 
equilibrium correlation function does not exist Q . How- 
ever, in the case of ordinary statistical mechanics Eq. (jSJ 



becomes equivalent to the Green-Kubo property. In this 
generalized sense we can state that the authors of Refs. 
E H II Q, H H studied the Green-Kubo relation of Eq. 
(O , all of them but the authors of Ref . [J , devoting their 
attention to the subdiffusional case. 

Herein we focus our attention on the Onsager relation 
of Eq. 10} . The only earlier work on the issue of the 
breakdown of the Onsager principle, caused by anoma- 
lous statistics, known to us, is that of Ref. 11]. However, 
here we plan to address the problem with the adoption 
of a Liouville-like approach, a fact that will allow us to 
establish some general conclusions concerning the break- 
down of the Onsager prescription. We shall illustrate 
the rules for the calculation of the four-time correlation, 
with a prescription that can be easily extended to cor- 
relation functions of any order. We expect that these 
prescriptions might lead to a successful evaluation of the 
fourth-order correlation function, which has been stud- 
ied so far by means of a factorization assumption which 
is violated by the non-Poisson statistics. 



II. AN IDEALIZED MODEL OF 

INTERMITTENT RANDOMNESS AND THE 

CORRESPONDING DENSITY EQUATION 

As done in Ref. U^, let us focus on the following 
dynamical system. Let us consider a variable y moving 
within the interval / = [0, 2]. The interval is defined over 
an overdamped potential V, with a cusp-like minimum 
located at y = 1. If the initial condition of the parti- 
cle is y(0) > I, the particle moves from the right to the 
left towards the potential minimum. If the initial condi- 
tion is 2/(0) < 1, then the motion of the particle towards 
the potential minimum takes place from the left to the 
right. When the particle reaches the potential bottom it 
is injected to an initial condition, different from y = 1, 
chosen in a random manner. We thus realize a mixture 
of randomness and slow deterministic dynamics. The 
left and right portions of the potential V{y) correspond 
to the laminar regions of turbulent dynamics, while ran- 
domness is concentrated at y = 1. In other words, this is 
an idealization of the map used by Zumofen and Klafter 
jl3!| , which does not affect the long-time dynamics of the 
process, yielding only the benefit of a clear distinction 
between random and deterministic dynamics. Note that 
the waiting time distribution in the two laminar phases 
of the reduced form has the same time asymptotic form 
as 



i;{t) = i^i-l) 



rpil- 



[t + TY 



(9) 



We select this form as the simplest possible way to ensure 
the normalization condition 



dtil){t) = I. 



(10) 



We note that Eq. (|10|l implies /i > 1. The condition 
fj, > 2 corresponds to the existence of a finite mean so- 
journ time, and, thus, to the possibihty itself of defining 
the stationary correlation function of the fluctuation ^, 
which, with the choice of Eq. ^ reads [lj| 



^dt) 






T 



t + T 



(11) 



From within the perspective of a single trajectory this 
dynamical model reads 



y = \[e{l-y)y^-Qiy-l){2~yr] + 



Ay{t) 



' random 



■Siy-i). 
(12) 

The function 0(x) is the ordinary Heaviside step func- 
tion, Aj, (i) is a random function of time that can achieve 
any value on the interval [—1,-1-1], and Trandom is the 
injection time that must fulfill the condition of being in- 
finitely smaller than the time of sojourn in one of the 
laminar phases. Note that z is a real number fitting the 
condition z > 1. In fact, the equality 



M 



(m-1) 



(13) 



relates the dynamics of Eq. H12() to the distribution 
Eq. ||UJ). The Poisson condition is recovered in the limit 
z — > 1, namely, in the limit /i ^ oo. Thus, in a sense, 
the whole region z > 1 {fj, < oo) corresponds to anoma- 
lous statistical mechanics. However, the deviation from 
normal statistical mechanics is especially evident when 
z > 1.5, a condition implying that the second moment 
diverges. In the case z > 2 the departure from ordinary 
statistical mechanics becomes even more dramatic, due 
to the fact that the first moment also diverges and, as we 
shall see in this Section, the process becomes non-ergodic. 
Let us move now to the density picture, namely, to a 
formulation of Eq. I|12|l from within the Gibbs perspec- 
tive. The form of this equation is: 



d 



piy,t) 



dt 
where 



-A|-[e(l-2/)y^ 
dy 



.e(y_l)(2-y)Xy,i)+C(i), 
(14) 



C{t) = 2\p{l,t). 



(15) 



It is important to stress that we are forced to set the 
equality of Eq. (|15|l to fulfill the following physical con- 
ditions 



dt ■ll=[0,2 



Piy,t)dy 



/=[0,2] 



^^piy,t)dy = 0, (16) 



which, in turn, ensures the conservation of probability. 
We assume the ordinary normalization condition 



7=[0,2] 



p{y,t)dy = 1, 



(17) 



which is kept constant in time, as a consequence of Eq. 
(|16|l . It is evident that the inhomogeneous term C{t) 
corresponds to the action of the stochastic term, namely, 
the second term on the right hand side of Eq. H12|l . 

It is important to point out that our dynamic per- 
spective allowed us to describe the intermittent process 
through the Liouville-like equation 



dpi^,y,t) 

dt 



T^p{^,y,t), 



(18) 



where y denotes a continuous variable moving either in 
the right or in the left laminar region, with ^ getting the 
values W or —W, correspondingly, and the operator TZ 
reading 



7^--A^[e(i-y)jy^ 



-eiy-l)i2-yr]+2\ / dyS{y~l). 
Jo 

(19) 

This operator departs from the conventional form of a 
differential operator, since the last term corresponds to 
the unusual role of an injection process, which is random 
rather than being deterministic. In the ordinary Fokker- 
Planck approach the role of the stochastic force is played 
by a second-order derivative, which is not as unusual as 
the integral operator of Eq. H19|l . Here the role of ran- 
domness is played by the back-injection process, which, 
from within the density perspective is described by an 
operator that selects from all possible values p{y,t), the 
specific value oi p{y,t) at y = 1. The idealization that 
we have adopted, of reducing the size of the chaotic re- 
gion to zero, with the choice of the process of back in- 
jection located a,t y — 1, makes it possible for us to use 
the continuous time representation and the equation of 
motion Eq. H18(l rather than the conventional Frobenius- 
Perron representation. This representation will allow us 
to obtain analytical results. However, the same physical 
conclusions would be reached, albeit with more extensive 
algebra, using the conventional maps and the Frobenius- 
Perron procedure described in the recent book by Driebe 

0. 

We note that the equilibrium probability density solv- 
ing Eq. ffT^ is given by 



Po{y) 



9(1 -y) , e(y-l) 



y 



i^-yy- 



(20) 



This equilibrium density becomes negative for z > 2 and 
signals the important fact that for z > 2 there no longer 
exists an invariant distribution. The lack of an invariant 
distribution accounts for the nonergodicity in the fluo- 
rescence of single nanocrystals, recently pointed out by 
Brokmann et al. [Ig. 

For the purposes of calculation in the next few Sec- 
tions, it is convenient to split the density p{y,t) into a 
symmetric and an anti-symmetric part with respect to 
2/ = l, 



p{y,i) ^ps{y,t) +PA{y,t). 



(21) 



This separation based on symmetry yields the following 
two equations from Eq. H14|) 

glPsiy.t) = -Xeil -y)-g-[y'ps{y,t)] 

+Xe(y-l)^[(2~yrps{y,t)]+C{t) 

(22) 
and 

-QlPAiy.t) = -Ae(l-2/) — [y>^(y,t)] 

+\Q{y~l)^[{2~yYpA{y,t)].{2i) 
dy 

We note that the anti-symmetric part of the density is 
driven by a conventional differential operator, which we 
denote by F. Thus, we rewrite Eq. (|23|l as follows 



d_ 
dt 



PAiy,t) ^tpA{y,t), 



where 



-Ae(l 



y)-y^-+Xeiy-l)-i2-yy 



(24) 



(25) 



The operator with the unusual form, containing C'(t), is 
only responsible for the time evolution of the symmetric 
part of the probability density. We notice, on the other 
hand, that any physical effect producing a departure of 
C{t) from its equilibrium value, if this exists, namely, if 
z < 2, implies a departure from equilibrium. A station- 
ary correlation function can be evaluated, as we shall see 
in the next few Sections, using only Eq. (|23|l . without 
forcing Eq. H22|l to depart from the equilibrium condi- 
tion. 

As we shall see in Section IVil the evaluation of corre- 
lation functions of order higher than the second cannot 
be done without producing a departure of C{t) from its 
equilibrium value. This might generate the impression 
that the correlation functions of order higher than the 
second cannot be evaluated without internal inconsisten- 
cies, if we use only the density picture. The evaluation of 
these higher-order correlation functions was done in Ref. 
[l7|. by using a procedure based on the time evolution 
of single trajectories. Actually, as we shall see in Section 
IVIl the density approach should yield the same result. 
However, we think that deriving this result using only 
the Liouville-like equation of this section is a hard task, 
which was bypassed in the past by means of the factor- 
ization approximation 12], violated by the non-Poisson 
case. 



III. THE CORRELATION FUNCTION OF THE 

DICHOTOMOUS FLUCTUATION FROM THE 

TRAJECTORY PICTURE 

Let us focus our attention on Eq. (|12|l and consider the 
initial condition j/o G [0,1]- Then, it is straightforward 



to prove that the solution, for y < 1 is, 

y{t) = yo{l~X{z-l)yr'ty'^'-^~'^ 



(26) 



From H2fci|) and imposing the condition y{T) — 1, we can 
find the time at which the trajectory reaches the point 
y = I, which is 



T = T{yo) = 



1 - yo~ 



Mz ~ i)yo" 



(27) 



Since we have to find ^(i) and ^(t) ~ ^{y{t)), from the 
general form of Eq. H27|l we obtain: 

^ = mi~yo)eiTiyo)-t) 

-6(2/0 -l)e(r(2-yo)-t)] 

+ 00 



Y^ sign 



i=0 




'i+i 



Q\Y.-^-t]-Q\Y,r,^t 



\k=0 



\k=0 



(28) 



where the time increments are given by 

ro - T(j/o)--i^^f4T©(l-yo) 

Mz-l)yo 



Ti>l = 






1-[1-A,(r0]- 
A(z-1)[1-A,(t,) 



-M^vin)). 



(29) 



Then, for the autocorrelation function we obtain the fol- 
lowing expression: 



{mm) 

VF2 



= ([e(l-yo)e(T(yo)-t) 

+e(yo-l)e(T(2-2/o)-0]) + 



+ 00 



^ ( sign(yo - l)sign 



i=0 



A. E 



Tfe 



^A;=0 



Q\Y.^k-t]-Q[Y.-^-t 



\.k=0 



\.k=0 



(30) 



As pointed out in Section ^I the calculation of the cor- 
relation function rests on averaging on the invariant dis- 
tribution given by Eq. (|20|l . As a consequence of this 
averaging, the second term in (|30|l vanishes. In fact, the 
quantity to average is anti-symmetric, whereas the sta- 
tistical weight is symmetric. 



It is possible to write the surviving term in the auto- 
correlation as 



(mm) 



= (2 - z) / e 



l-y'- 



1 



'^\X{z-l)y^-^ Vy-i"^^ 

n(l+A(2-l)t)-l/("-l) 

= (2 - z) / y-'+'dy 



with 



= (l + A(z-l)t)-(2--)/(--i) 



P = 



(32) 



Since we focus our attention on < /? < 1, we have to 
consider 3/2 < z < 2. Note that the region 1 < z < 3/2 
does not produce evident signs of deviation from ordinary 
statistics. However, as we shall see in Section|V] an exact 
agreement between density and trajectory is recovered 
only at z = 1, when the correlation function becomes 
identical to the exponential function exp(— Ai). Note also 
that Eq. H31|l becomes identical to Eq. (|ll|l after setting 
the condition 



A(z-l) 



1 
T' 



(33) 



IV. THE CORRELATION FUNCTION OF THE 

DICHOTOMOUS FLUCTUATION FROM THE 

DENSITY PICTURE 

The result of the preceding Section is reassuring, since 
it establishes that the intermittent model we are using 
generates the wanted inverse power law form for the cor- 
relation function of the dichotomous variable ^(i). In this 
Section we show that exactly the same result can be de- 
rived from the adoption of the Frobenius-Perron form of 
Eq. (O. 

To fix ideas, let us consider the following system: a 
particle in the interval [0, 1] moves towards y = 1 follow- 
ing the prescription y = Xy^ and when it reaches y — 1 
it is injected backwards at a random position in the in- 
terval. The evolution equation obeyed by the densities 
defined on this interval is the same as Eq. (|14|l . with 
C{t) ~ Xp{l,t). This dynamic problem was already ad- 
dressed in Refs. [l8lll9|. and solved using the method of 
characteristics as detailed in Ref. |2fl] . It is important to 
stress that this approach is the requisite price for adopt- 
ing the idealized version of intermittency. The adoption 
of the more conventional reduced map of Ref. [l3| would 
have made it possible for us to adopt the elegant pre- 
scriptions of Driebe [T^ , as done in Ref. [13 ■ 

Let us remind the reader that the solution afforded by 
the method of characteristics, in the case of this simple 
nonlinear equation with stochastic boundary conditions 



P{v,i) 



Ap(l,0 



—,d^ 



9M-^)y'-^) 
y 



tp 



[9Ay'-'t)V/^ 



X 



1 



9ziy'-'t) 



(31) where 



5,(a:) = (l + A(z-l)a;)^/(^-i). 



(34) 



(35) 



To find the autocorrelation function (C(i)C(O)) using 
only densities, we have to solve Eqs. H22() and H23() . which 
are equations of the same form as that yielding Eq. H34|l . 
For this reason we adopt the method of characteristics 
again. To do the calculation in this case, it is convenient 
to adopt a frame symmetric with respect to y = 1. Then, 
let us define Y — y — 1. Using the new variable and Eq. 
(IMJI . we find for Eqs. if^ and if^ the following solution 



Ps{Y,t) 



Aps(0,r) 



g,{{t-T){l-\Y\y-^) 

i-iri 



dT 



PS 



1 



[gA{l-\Y\Y-H)]y^ 



1 



gA{l-\Y\r-H) 



,0) X 
(36) 



and 

PA{Y,t) - p^(sign(y) 
1 



1- 



l-ll^l 



M{i-\Y\r-H)]^/^ 



g{{l~\Y\r-Hy 



,0 1 X 

(37) 



Then, the solution consists of two terms: (1) the former 
is an even term and is responsible for the long time limit 
of the distribution evolution and (2) the latter is an odd 
term which disappears in the long-time limit. We note 
that (2) is a desirable property because Eq. H23|l does 
not contain the injection term C'(t) and the equilibrium 
density H2U|I is an even function, independently of the 
symmetry of the initial distribution. 

As pointed out in Section^ the two-time correlation 
function is determined by the anti-symmetric part of the 
probability density alone. Thus, the definition of the 
autocorrelation function does not conflict with the equi- 
librium assumption. In Section we shall see that this 
conflict reemerges when we attempt to evaluate higher- 
order correlation functions. Let us calculate the autocor- 
relation function explicitly: 



(mm) = (2-^) 



1 



[i-ry 
1 



T = l- 



la^{(i-y)'~^t)]^^' 



■g.ia-Yy-H) 



dY 



(38) 



The integral H38() is exactly solvable and leads to the ex- 
pression 

(mm) - (1 + A(z - i)t)-(2-^)/(^-i', (39) 

which is the same result as that found in Section IIIII 
using trajectories rather than the probability densities. 

In a similar way, it is possible to calculate the auto- 
correlation {Y{t)Y{Q)) and determine that its temporal 
behavior is an inverse power law with the same exponent 
as that in Eq. l|5!?|) . 



V. ONSAGER REGRESSION TO EQUILIBRIUM 

In conclusion, in the two preceding Sections we have 
established that the Liouville-like representation of Eq. 
(|18|l yields, as expected, the equilibrium correlation func- 
tion of Eq. lfTT|l with 



T = 



A 



(40) 



What about Onsager's regression to equilibrium? Let 
us first of all discuss a physical condition where the pre- 
scription of Eq. Q) is fulfilled. Let us consider the initial 
distribution p(y, 0) defined as follows 



and 



p(y,0) = 0, y<l 



piy,0) = 2poiy), y > 1, 



(41) 



(42) 



where po{y) denotes the equilibrium distribution of Eq. 
(|20|l . It is convenient to point that this equilibrium dis- 
tribution is symmetric and that the factor of 2 serves the 
purpose of normalizing the out-of equilibrium condition 
that we are studying. Let us denote by AP{t) the popu- 
lation difference between the left and the right state. The 
choice we made sets the initial condition AP(0) = 1. It is 
straightforward to show that Onsager's regression implies 
that AP(t) = <&5(t). Let us split the initial distribution 
in the symmetric and anti-symmetric part. 



P(y,0) =pa{y) 



-P^r\y,o)- 



We note that 



and 



p'x'\y,o)^~Poiy),y<i 



p'i'\y,0)^Po{y),y>l. 



(43) 



(44) 



(45) 



The symmetric part does not contribute to AP{t), only 
the anti-symmetric part does. Thus, we obtain 



AP{t) = 2 



/ cxp(fi)p' 
Jo 



^r^iy.o), (46) 



where F is the operator driving the anti-symmetric part, 
defined by Eq. H25|) . With the choice of initial condition 
made it is straightforward to prove that the right end 
side of Eq. (|46|l is the equilibrium correlation function 
evaluated in Sections IIIII and Hvl Thus, the prescription 
of ordinary statistical physics is fulfilled. 

What about the regression to equilibrium in general? 
We note that we can adapt the earlier arguments to any 
initial condition thereby yielding 



AP{t) 



exp(R)p^"°"^')(y,0). 



(47) 



Note that p^°"'^'^ {y,t) obeys the Liouville-like prescrip- 
tion corresponding to -jI = Ay^, namely, the equation 
of motion of the left laminar region, without any back- 
injection process. Since, as we have seen in Section 
im the back-injection term serves the purpose of keep- 
ing constant the population of the system, we have that 
AP(t) -^ for t -^ oo. More precisely, we obtain 



dAPit) 
Jt 



2r / exp(ff)p^-"^«)(y,0) 



-2p^"°"^'^(l,i). 

. ^^^^ 
The superscript noneq serves the purpose of pointing 

out that in general the perturbation process creating the 
necessary initial asymmetry will not realize, for the left 
portion of the asymmetric component, a form exactly 
identical to the left portion of the equilibrium distribu- 
tion. This observation makes it possible to estimate the 
time asymptotic behavior of AP{t) in general. The exact 
time evolution of AP{t) depends on the detailed effects 
of the perturbation. However, for any realistic pertur- 
bation, we can prove that the time asymptotic behavior 
obeys an universal prescription. If the perturbation does 
not affect the equilibrium distribution in the regions close 
to the borders, namely for y < e and y > 2 — e, where 
e ^ 1, we have that p^^ {y-, 0) vanishes for y < e. To 
evaluate the asymptotic behavior of AP{t) we select the 
infinitesimal portion dy{0) of the interval [0,1], closest 
to J/ = 0, where p^^""'^*^ (y^ 0) does not vanish. We call 
M the number of trajectories located in this interval at 
t = 0. The asymptotic behavior of AP{t) is determined 
by the time necessary for these trajectories to reach the 
border y = 1. The first trajectory will reach the border 
after a given time t — T first, after which the number M 
will begin decreasing, thereby determining the decay of 
AP{t) in this time asymptotic region. It is straightfor- 
ward to prove that the time of arrival at y = 1 of the 
trajectory with initial condition y(0), called t, is related 
to y(0) by 



2/(0) = ^ 

[l + (z-l)Ai]- 

Thus we obtain 

dAP{t) dM 1 



(49) 



dt dt [1 + (z - l)\t] 



^- (50) 



By integrating this equation and taking into account that 
z = ^/{^ — 1), we finaUy obtain 

hm AP(t) ex -1^, (51) 

which sanctions the breakdown of Eq. Q . 



VI. HIGHER-ORDER CORRELATION 
FUNCTIONS 

We now show that the calculation of higher-order cor- 
relation functions, though difficult, can be done by es- 
tablishing an even deeper connection with the quan- 
tum mechanical perspective. Let us address the prob- 
lem of evaluating the fourth-order correlation function 
(C(*4)C(^3)C(^2)C(^i))- According to the prescription that 
we adopted in Section llVl we must proceed as follows. We 
move from the equilibrium distribution and let it evolve 
for a time ti. The distribution selected is in equilibrium. 
Therefore it will remain unchanged. At time i = ti we 
apply the operator ^ to the distribution. Since the equi- 
librium distribution is symmetric, the application of the 
sign operator changes it into the anti-symmetric distribu- 
tion. We let this distribution evolve for the time ^2 —ti- 
This has the effect of yielding 



^dt2-h)p'^r\y^^) 



(noneq) / , 

Pa (y>*2 



h). 



At this stage there are two possibilities: 



S' 



eq) 



(a)pr^^^^(y,i2-ti) = 0; 

(b)pi"°"^')(y,i2-ti)^0. 

Let us consider the case (a) first. In this case, we pro- 
ceed as follows. At time t2 we apply to the distribution 
the operator ^, and we change it into the original equi- 
librium distribution. This means that the time evolution 
from ^2 to ty, leaves it unchanged. At time t^ we apply to 
it the operator ^ and we turn it into p^'^ again. We let 
this distribution evolve till to time ^4. At this time we 
apply to it the operator ^ again and we make the final 
average. The result of the condition (a) yields 

(C(t4)c(t3)c(t2)c(ti)) = {^{unh))m2nh)). (52) 

In a recent work 17] it has been shown that this fac- 
torization condition is violated by non-Poisson statistics. 
The demonstration was made by applying the method of 
conditional probabilities to the study of single trajecto- 
ries. Thus we are forced to consider case (b). 

The problem with condition (b) is that it yields a dis- 
tribution with the symmetric component departing from 
equilibrium. This departure is in an apparent conflict 
with the assumption that the autocorrelation function 
is calculated using the equilibrium condition. In fact it 
seems to be equivalent to stating that the calculation of 



an equilibrium correlation function generates an out of 
equilibrium condition. Let us see why. 

In the case of the two-time correlation function we ap- 
ply the operator ^ to the density distribution twice. The 
first application allows us to observe the time evolution of 
the anti-symmetric component distribution density, with 
no conflict with equilibrium, given the fact that the sym- 
metric component remains at equilibrium. The second 
application of ^ turns 



$5(^2 



into 



ti)p^r\y,Q) 



(eg) 



(noneq) 

■Pa 



iy,t2 -h) 



^dt2-h)pr'{y,0) 



(noneq) / , 

Ps (y>*2 



ti). 



The calculation done in Section IIVI proves that 

T^T^lPg^iy, t2 — ti)] = 0, with the symbol Tr denoting, 
for simplicity, the integration over y from to 2. To eval- 
uate the fourth-order correlation function, after applying 
the operator S, for the second time, we must study the 
time evolution of Pg°"'^*(y,t2 — ii) from ^2 to is, yielding 

to pg"^'^'{y,t3,t2,t2 — ti). This is a contribution gen- 
erating some concern, since it activates again the back 
injection process, which we have seen to be intimately 
related to the deviation from equilibrium. However, the 
compatibility with equilibrium condition is ensured by 

the property Tr[pl?°""''^(y, ^3,^2,^2 - ^i)] = 0. 

At time ^3 we have to apply the operator ^ again, and 
this allows us to make an excursion in the anti-symmetric 
representation, with the time evolution given by the op- 
erator exp(r(i4 — t^)). At time ti we apply the operator 
^ again, we go back to the symmetric representation and 
we conclude the calculation by means of the trace oper- 
ation. 

In conclusion, the compatibility with equilibrium is 
guaranteed by the fact that at the intermediate steps of 
the calculation Tr[p] = (the final step, of course, gener- 
ates the density-generated correlation function, thereby 
implying Tr[p] ^0). If the intermediate p is anti- 
symmetric, this condition is obvious. If the intermediate 
p is symmetric, the vanishing trace condition generates 
the apparently unphysical property that the symmetric 
contribution gets negative values over some portions of 
the interval /. We have to stress, however, that the 
Liouville-like approach illustrated in this paper keeps the 
distribution density p{y, t) definite positive, as it must. 
The generation of a negative distribution density refers to 
the calculation of the equilibrium correlation functions, 
of any order, and it is a quantum-like property that must 
be adopted to guarantee that the genuine distribution 
density remains in the equilibrium condition. 



VII. CONCLUDING REMARKS 

In conclusion, this paper shows how to derive the equi- 
librium correlation function using only information af- 
forded by the Liouville-like approach. The intriguing 
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problem to solve was how to use the Liouville equation, 
without conflicting with the equilibrium condition. The 
solution of this intriguing problem is obtained by split- 
ting the Liouville equation into two independent com- 
ponents, the symmetric component corresponding to Eq. 
()22|l and the anti-symmetric component corresponding to 
Eq. H23|) . This splitting allows us to study the regression 
to equilibrium of the correlation function $ ^ (t) , through 
Eq. H23|l . without ever departing from equilibrium, a 
physical condition that is controlled by the independent 
equation of motion for the symmetric component, Eq. 
(|22|l . This is very formal, and Section^makes it possible 
for us to reconcile it with physical intuition. We imagine 
that the Liouville equation is used to evaluate the differ- 
ence between the population of the left and right state. 
This makes it possible to establish a direct connection 
between the experiment of regression to equilibrium and 
the formalism of SectionslTlland irvl We have to create an 
asymmetrical initial condition, with more population on 
the left than on the right. The time evolution of AP(i) 
depends only on the time evolution of the anti-symmetric 
component of the distribution density, and consequently 
only on the operator F. This is the reason why it is pos- 
sible in principle to connect regression from an out of 
equilibrium initial condition to the equilibrium correla- 
tion function, a condition that implies no deviation from 
equilibrium distribution. However, creating in a finite 
time an out of equilibrium condition such that the left 
part of the anti-symmetric component is identical to the 
left part of the equilibrium distribution, is impossible. 
This is the reason why we predict the breakdown of the 
Onsager principle in general. 

Section lvTI cxplains why in the literature on dichotomic 
fluctuations the factorization assumption of Eq. I|52|) is 
often made regardless of the Poisson or non-Poisson na- 
ture of the underlying process. See, for instance, the work 
of Fulinski |2j] as well as Ref. [l3|- In fact, if condition 
(a) applies, the higher-order correlation functions are fac- 
torized, thereby making their calculation easy. However, 



this assumption conflicts with the trajectory arguments 
of Ref. ,131 J which prove that the factorization condition 
is violated by the non-Poisson condition. A rigorous use 
of the Liouville equation shows that condition (a) does 
not apply, and that we have to use condition (b) instead. 
The calculation of the fourth-order correlation function 
is not straightforward, and this is the reason why, to the 
best of our knowledge, it was never done using the Liou- 
ville approach. 

This is a fact of some relevance for the creation of mas- 
ter equations with memory. There are two major classes 
of generalized master equations. The first class is dis- 
cussed, for instance, in Ref. [Tl|. The master equations 
of this class are equivalent to the Continuous Time Ran- 
dom Walk (CTRW) of MontroU and Weiss 1\ and are 
based on the waiting time distribution ■0(t). The sec- 
ond class of master equations is based on the correlation 
function $^ (t) . Recent examples of this second class can 
be found in Ref. [I^l and in Ref. 22J. Due to the direct 
dependence on the correlation function ^^{t), the deriva- 
tion of the master equations of this second class is made 
easy by the factorization assumption. It must be pointed 
out, on the other hand, that the factorization property, 
which is not legitimate with renewal non-Poisson pro- 
cesses, is a correct property if the deviation from the 
exponential relaxation is obtained by time modulation of 
a Poisson process 221 ■ Beck 23j is the advocate of the 
modulation process as generator of complexity. Thus, we 
find that the master equations of the first class are gen- 
erated by the renewal perspective of MontroU and Weiss 
and those of the second class are the appropriate tool 
to study complexity along the lines advocated by Beck 
[23 ■ We think that the results of the present paper might 
help the investigators in the field of complexity to make 
the proper choice, either modulation or renewal 24P|, or a 
mixture of the two conditions. 
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